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Abstract Consider a multidimensional SDE of the form Xt = x + J* b(X s -)ds + 

J'f(X,-)dZ. where (Z s ) s >0 is a symmetric stable process. Under suitable assump- 
tions on the coefficients the unique strong solution of the above equation admits a 
density w.r.t. the Lebesgue measure and so does its Euler scheme. Using a parametrix 
approach, we derive an error expansion w.r.t. the time step for the difference of these 
densities. 
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1 Introduction 

Consider the following ]R d -valued Stochastic Differential Equation (SDE in short) 

X t =x+ [ b(X s -)ds+ [ f(X s -)dZ s , (1.1) 
Jo Jo 

where b, f are respectively Lipschitz continuous mappings from R d to M. d and M. d to 
M. d Cg) R d and (Z s ) s >q is a general Levy process. The previous assumptions guarantee 
the existence of a unique strong solution to (|l.ip . Also, this solution satisfies the strong 
Markov property, see e.g. Theorem 7 and 32 Chapter 5 in Protter |Pro04| . Let T > 
be a fixed time horizon and {X^) t ^j\ a given approximation scheme of (-Xt)tg[o,Tl 
associated to the time step h := T/N, N G N* on the grid A := {tj := ih, i G 
[0, NJ}. When speaking about weak approximation of {O} two kinds of quantities are 
of interest. The first one writes 

£l(x,T,N) ~E x [g(X T )]-E x [g(x¥)] 
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for a suitable class of test functions g. The second one concerns, when it exists, the 
approximation of the transition density p of the original SDE If the approximation 

scheme (X^) tej \ admits as well a transition density p N , the quantity under study 
becomes 

£ 2 (x,y,T,N) : = (p — p N )(T,x,y). 

In both cases, the goal is to give a bound or an error expansion of these quantities 
in terms of h. The error expansions are particularly useful for practical simulation. 
For £ i , the expansion allows to use the Romberg Richardson extrapolation to improve 
the convergence of the discretization error, see e.g. Talay and Tubaro |TT90| . On the 
other hand, if p and p N exist, and a suitable expansion of £2 holds, it can be useful to 
estimate the sensitivity of £\ w.r.t. to the spatial variable x and it also allows to get a 
control on £\ for a wider class of test functions g than those considered by the direct 
methods used to control this quantity, see e.g. Guyon Gu y06| . Indeed, the typical 
assumptions and techniques associated to the study of £ 1 and £2 are quite of different 
nature. 

In the continuous case, i.e. Z s = ba + aW s where {W s ) s >o i s a standard d- 
dimensional Brownian motion, provided the test function g and the coefficients b, f 
are sufficiently smooth and g has polynomial growth, without any additional assump- 
tion on the generator Talay and Tubaro |TT90| derive an error expansion at order 1 
for £i(x,T, N) when (X^) t&J \ is the Euler approximation. Their proof is based on 
standard stochastic analysis tools: Ito's expansions and stochastic flows. To obtain the 
same kind of result for bounded Borel functions g some non degeneracy has to be 
assumed, namely hypoellipticity of the underlying diffusion, and the proof relies on 
Malliavin calculus techniques, see Bally and Talay |BT96aj . The authors also manage 
to extend their results to £2(2, y, T, N) for a slightly modified Euler scheme |BT96bj . 

Anyhow, in the uniformly elliptic case, the most natural approach to handle the 
estimation of the quantity £2(x,y,T,N) consists in using the so called "parametrix" 
technique introduced to obtain existence and controls on the fundamental solutions of 
PDEs, see e.g. Mc Kean and Singer |MS67| or Friedman [Fri64| . Roughly speaking it 
consists in expressing the density of Xj- in terms of an infinite sum of suitable iterated 
kernels applied to the density of an SDE with constant coefficients. This has been done 
successfully by Konakov and Mammen |KM02j .The main advantage of this approach is 
that the density of the solution Xj- and the Euler approximation X^ can be expressed 
in the same form and therefore quite directly compared. Furthermore this technique 
turns out to be quite robust and can be applied as soon as good controls on the densities 
p,p N and their derivatives are available, sec e.g. KMM08 for an extension to a slightly 
degenerate framework. 

For a general Levy process Z and suitable smooth functions b, f, g, under additional 
assumptions on the behavior at infinity of the Levy measure v of Z, that is integrability 
conditions of the large jumps, Protter and Talay |PT97j . manage to get a control at 
order one or even an error expansion for £\(x,T, N) with the same approach as in 
|TT90| . In that work the approximation is the Euler scheme which for a general Levy 
measure v cannot always be exactly simulated on a computer. 

The quantity £ \(x, T, N) for approximations of the Euler scheme that can be sim- 
ulated has also been studied by Jacod et al. [JKMP05] who derived bounds at order 1. 
Moment conditions are also assumed. We finally refer to the work of Hausenblas and 
Marchis |HM06] for approximations of Poisson jump measures that are easy to simulate. 
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In this work, we consider the case where (Zt)t>o i s an ^-stable symmetric process, 
a G (0, 2). Under suitable non degeneracy assumptions on its coefficients specified below 
(see (A-l)-(A-3)), equation 1)1. 1|) is known to have a density p w.r.t. the Lebesgue 
measure. This can be proved via a Malliavin calculus-Bismut integration by parts 
approach, see e.g. Bichteller el al. |BGJ87| . Also, a direct construction of this density 
using a parametrix expansion has been obtained by Kolokoltsov [KolOO] who derived 
as well "Aronson's like" bounds with time singularity depending on the index a of the 
stable process (Zt)t>o- 

Analogoulsy to the "diffusion case" the first step of the parametrix is to consider 
that the density p(T, x, y) of JXTTJ) can be approximated by the density of the process 
Xf = x + b(y)t + f(y)Zt at time T. Namely, we freeze the coefficients in at the 

final spatial point. The next crucial point is to obtain sharp estimates of the stable 
density p^{T,x, .) of Xj, and its derivatives in order to solve the parametrix integral 
equations. 

Stable driven SDEs appear in various applicative fields, from mathematical physics 
to electrical engineering or financial mathematics, see IP06 , SK74 or [JMW96], there- 
fore their approximation becomes of interest. To approximate equation Ql.ip . setting 
<j>(t) := inf{£; : t j < t < f;+i}, we introduce the Euler scheme 

X? = x + [ b(X% (s) )ds+ [ f(X^ s) )dZ s . (1.2) 
Jo Jo 

The computation of the above scheme only requires to be able to simulate exactly the 
increments of (Zt)t>o> which up to a self similarity argument only amounts to simulate 
a stable law. This aspect is for instance discussed in Samorodnitsky and Taqqu [5T94 , 
Weron and Weron |WR95| or Section 3 of [PT971 - Under the same assumptions (A- 
l)-(A-3), the Euler scheme defined in Q1.20 also has a density p N . 

Observe that the results of |PT97] . [JKMP05] cannot be directly applied, even for 
the study of £\ (x, T, N), since stable laws have heavy tails. Comparing the parametrix 
developments of p and p N we obtain an expansion with leading term of order 1 in h for 
£2(2;, y, T, N). The parametrix expansion of p is discussed in |Kol00) . see also Section 
and Appendix, whereas the parametrix expansion of p N can be related to the ideas 
developed in |KM00 , K M02] for the diffusive case corresponding to an index of stability 
equal to 2. 

This result also emphasizes the robustness of the method that naturally extends to 
a broad class of processes. Let us mention that, using a Malliavin calculus approach, 
Hausenblas |Hau02] , derived an upper bound of order one w.r.t. h for the quantity 
S\(x,T, N), g £ L°° in the scalar case. Concerning functional limit theorems for the 
approximation of stable driven SDEs we refer to the work of Jacod |Jac04| . 

The paper is organized as follows. In Section [2] we state our standing assumptions 
and main results. In Section we prove the existence of the densities for both the 
stable driven equation and its Euler scheme and also give a parametrix representation 
of these densities. Section [4] is dedicated to the proof of the main results. Eventually, 
we state in Section [5] weaker assumptions under which our main result holds and we 
also briefly discuss how to extend it to the case of a stable process perturbed by a 
compound Poisson process. 
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2 Assumptions and Main results 

2.1 Assumptions and Notation 

In the following we consider symmetric stable processes, that is, for all t > 0, u G K d , 

E[exp(z( W , Z t ))] = exp(ii( 7 , »)+tf ~ 1 - ^fe^*")) 

= exp(tt( 7l «>-t / |<«,u)| a A(ds)), (2.1) 

where A is a symmetric measure on the unit sphere 5 d_1 (i.e. for every A in the 
Borel cr-field B(S d ~ 1 ), X(A) = \(-A)). The second equality in equation (|2~Tj) is then 
obtained by direct integration over p and A = C Q A with 

, . _i f Tva\ TV 

C a --r(l-a)a cos I — I l a -tx + 2 1I "=i- 

We refer to the proof of Theorem 9.32 in Breiman |Bre68| and Lemma 2, Chapter 
XVII. 4 in Feller |Fel66j for the expression of C a . 

We now introduce our assumptions. Fix an integer q > 2. We assume that 

(A-l) For d > 2, the spherical measure A has a C q (S d ~ 1 ) surface density and for all 
d > 1, there exist constants < C\ < C2 < +00, Vp G R d , 

d\p\ a < f \(p,s)\ a X(ds)<C 2 \p\ a . 
Js d - 1 

(A-2) The coefficients b and / and their derivatives up to order q are uniformly 
bounded in x. Thus, for 1 < a < 2, B(x) := b(x) + f(x)j is uniformly bounded. 
We impose for < a < 1, B(x) = for all x G R d . 

(A-3) There exist constants < c < c < +00 s.t. for all x G R d , £ G R d , 

ci£i 2 <</(z)e,o<ci£i 2 . 

From now on we assume that Assumptions (A-l)-(A-3) are in force. 

Remark 21 Note that for d = 1, with the convention S° = { — 1, 1}, we have C\ = C2 
in (A-l) even without symmetry. The symmetry is actually not needed in that case, 
see the beginning of Section 3 in JKolOff/ . 

Remark 22 The zero drift condition in (A-2) comes from the fact that for a G (0, 1] 
the addition of a drift of order t does not correspond to a negligible term in small time 
with respect to the natural scale t 1 ^", see Appendix B in \KM1Q$ for details. 

In the following we denote by C a positive generic constant that can depend on a, d, 
the bounds appearing in the previous assumptions but neither on nor on the spatial 
points involved. Its value may change from line to line. Other possible dependencies, 
especially w.r.t. the final time T are explicitely specified. Concerning functional spaces, 
we denote by C^(R d ), k G N*, the Banach space of continuous bounded functions 
having bounded derivatives up to and including the order k with the norm ||/|| := 
m&x 0<l<k sup xm d |/")(a;)|. Eventually Co(K d ) stands for the functions in C^(R d ) 
with compact support. 
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2.2 Generator 

From equation (|2.ip and standard computations, see e.g. equation (5.11) in JKMP05J, 
we derive that for every smooth function g G Cg(M d ), the generator of (|1.1[) writes 

$g{x) = (B(x), V x g{x)) - f g(x + f(x)y) - g{x) - {Vx9{ f ^ )v) v(dy), 

where B(x) — b(x) + f{x)~y and v stands for the Levy measure of Z. Introduce for all 
A G B(R d ), Uf^^A) := u({y G R d : f(x)y G A}) and denote by XfM its spherical 
part (which is still a symmetric measure). Setting z — f(x)y in the above equation, 
using the symmetry and the polar coordinates we derive: 

&g(x) = {B(x),X7 x g(x)) + 
/ s _ jf" (*(» + /»)"*(») - ^A /( .)W- (2-2) 

Remark 23 Denote similarly to H2.1|) , Xf( x ) = C a Xf( x )- The uniform ellipticity con- 
dition (A-3) allows to have good controls on the measure \f/ x \(-). As a consequence 
of (A-l), (A-3) one gets that there exist constants < C_ 1 = Cj(c, d, ct) < C2 = 
C 2 (c, d, a) < +00 s.t. Vp G R d , x G R d , 

CM a <l \(p,s)\ a \ f{x) (ds)<C 2 \p\ a . (2.3) 

2.3 Main results 

Proposition 21 For every t > the solution Xt (resp. ) of (|1.1[) (resp. (|1.2|l ) 
has a density p(t,x,-) (resp. p N (t,x,-)) w.r.t. the Lebesgue measure. Additionally, as 
a function of the space variables the density p is in C*(R d x R d ) if a > 1 and in 
C^~ 1 {R d x R d ) ifa<\. 

To state the theorem we first need some notation. Introduce for all £ G R d and all 
smooth function ip(t, x, y) the integro-differential operators: 

$ ( <p(t, x, y) = V x <p{t, x, y)) + 

L C + "> - *> - "^rf^) ;&*/»«■>■ (i4) 

With this definition we write for given (x,y) G M. d : 

$*<p(t,x,y) = $ v <p(t,x,y), Mm GN*, h£*J p(t, s, y) = f * c j <p(t,x,y) |f= x , 

(2.5) 

Note that we have 4>*tp(t, x,y) = $ip(t,x,y) defined in (|2.2p but in general, for m > 2, 
<p{t,x,y)^{3>) m <p{t,x,y). 
Define now, for t > 0, the kernel 

H(t,x >V ):=(*-*„)p , '(t,x,») (2.6) 
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where p^it, x, y) denotes the density at point y of Xt = x + b(y)t+ f(y)Zf. Note that the 
variable y acts here twice: as the argument of the density and as a defining quantity of 
the process Xt(= Xt.x.y), i-e. the coefficients are frozen in y. Eventually we introduce 
the continuous and discrete convolution operators 



jf"*/ 



dzip(u,x,z)ij)(t — u, z,y),Vt £ [0, T], 



<p ®n *P(t,x,y) = I du I dzip(cj>(u),x,z)ip(t - cj>(u), z,y),Vt e {(^)ie[i,JV]}> 

with cf>{u) is defined just before (|1.2p and denotes the largest discretization time lower 
or equal to u. Also ip g! = </? and ip ® ff( r ' = ^ ® i/( r_1 ^ ® stands for the 
r-fold convolution. 

Theorem 21 Suppose q > c! + 4. Tafce < M < q — (d + 4). There exists a function 
Ritf(T,x,y) with \R M (T,x,y)\ < C M {T) y l+ \ y } x \ d+ S ) ~ Pa,Ai(T,y ~ x ) f° r some 
positive constant Cm(T) such that 



(p-p N )(T,x,y) = Yl 



1=1 



(i + iy. 



i+i 

p®N ) p 



(T,x,y) 



E 



fc=l 



(fc + 1)! 



P ®N 



fc+1 



iV 



(T, y) + h M R M (T, x, y) 



th Ej=i-- - = ^ Vt G {(*i) <6 [i,jv]}, P d (t,x,y) := J2T=o (p®JV^ W ) 



M-l 

E 

2=1 



(+1 



M-l 

E 



fc+i 



p a ®jv(^ -**] P W ) (T,x,y) 



< Pa,Ai(T,y - x), 

< Pa,Ai(T,y- x). 



Remark 24 In the above expression, one writes for all I £ [1, M — 1], 



(+1 



!+l 



$-$*) <p(t,x,y) = ^Cf +1 * fc (-^)' +1 -V(*,x,y), 



fc=i 



whereas, Vfc £ [1, M - 1] 
fc+l 



/ \ K-f-l 



times 

(* c -*„) fc+ V(*,a!,v) 
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Remark 25 The terms in the previous expansion depend on N . Anyhow using itera- 
tively the Theorem and controls on <8>at — ® (see also Lemma \4~l\ ) it is possible to obtain 
an expansion with terms independent of N . For small M explicit formulas are thus 
easily derived but in all generality the terms become less transparent. For M = 2 one 
gets 

(p - p N )(T, x,y) = \ (p ® N (* - $*) 2 p d - / ® N (*. - (T, x, y) 

h ( ~* o ~ ~* +Ji 2 R 2 (T,x,y) 

= g (P ® (* - * ) P - P ® (** " * ) Pj ( T > K > 2/) 

+h 2 R2(T,x,y) = ^(p® (<£ 2 -#2)p)(T,«,y) + /^(r, z, y), 

where R,2(T, x, y) < C(T)p a ^(T, y — x) for some positive constant C(T). 

From the above expansion and the controls on the density and its derivatives, see 
e.g. Theorems 3.1, 3.2 and Proposition 3.1 in fKolOOf or Lemma pg| we can derive 
the error expansion for £\{x,T, N) for measurable functions g satisfying the growth 
condition 3C > 0, \g(x)\ < (7(1 + |x|^),/3 < a. In particular, we do not need the 
smoothness assumption on g required in the approach of JTT90L !PTW^. We recall 
that the expansion of £\{x,T,N) allows from a practical point of view to improve the 
convergence rate of the discretization error using the Romberg Richardson extrapolation. 
This simply consists in observing that the expansion yields E[g(X^p)] — (K[2g(X^ N )] — 
K[g(Xj()]) = 0(h 2 ). The associated Monte Carlo estimator, involving a refined scheme, 
is then used for simulations see JTT9ff/ for details. 

Also, the expansion can be used to study the sensitivity of £i(T, x, N) w.r.t. x with- 
out any additional assumption on g. This is crucial for financial applications (hedging) , 
see e.g. Guyon \Guy06\j for further developments in the diffusive case. 



3 Stable driven equations and their Euler scheme: existence of the density 
and associated parametrix expansion 

3.1 Stable driven equation 

3.1.1 Proof of Provosition \21\ existence of the density for the solution of 1)1. 1|) 

For (Xt)t>o, the existence of the density derives from Proposition 3.4 in |Kol00| . where 
some properties of the fundamental solution of dtp(t,x,y) = <l>p(t,x,y),p(0,x,y) — 
5(y — x) are discussed, and a standard identification argument, see e.g. Dynkin |Dyn63| , 
Theorem 2.3, p. 56. The stated smoothness of the density is then a consequence of point 
(ii) of the same Proposition. 

Remark 31 The existence of the density is discussed in Bichteler et al. \BGJ81^ . 
where it is proved thanks to a Bismut-Malliavin approach. This technique requires the 
computation of a tangent equation associated to the gradient flow that involves the 
derivatives of the coefficients of equation {TTTJ . Thus, some additional smoothness of the 
coefficients is needed, see e.g. Theorem 6.48 of the above reference. We also mention the 
result of Picard \Pic9ty . Theorem 4-1, that gives existence and smoothness of the density 
for Levy driven SDEs for very singular Levy measures, provided there are sufficiently 
small jumps. For smooth coefficients b, f, it includes in particular the case of (|1.1|) 
where the spherical measure X can be atomic. 
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3.1.2 Parametrix expansion of the density 

For the sake of completeness and also because it is crucial for the discrete model we 
briefly recall how to get through a "parametrix" approach a series expansion for the 
density p(t, x, y). 

Introduce, for all x, y £ K d the following stochastic "frozen" stable driven equation 
Xt = Xt.x.y defined for t > by 

X t =x+ f b(y)du + f f(y)dZ u . (3.1) 
Jo Jo 

By computation of the Fourier transform of Zt and Fourier inversion the transition 
density jP(t, x, z) of Xt at point z £ M. d explicitly writes 

p»(i,x J *) = -i- f J e-^-'^^expj-t J<p )S )rA /fo) (d S )Jdp,(3.2) 

where ^f( y ) has been introduced in Section \2. 2 1 The densities of the solutions of (13.111 
and (|1.1[) satisfy respectively 

^-(t,x,z) = $ vP ^(t,x,z) y for t > 0, Or, 2) £ (R d ) 2 , p»(0, x, z) = S(z - x), 
^(t,x,z) = $p(t,x,z), for t > 0, (z,z) £ (R d ) 2 , p(0,x,z) = 5{z-x). (3.3) 

Note carefully that the derivatives in <& y are taken w.r.t. the x variable. 

We will speak about the operators appearing in (|3.3p as the "frozen" and "unfrozen" 
ones. In the following V(t, x, z) £ R + * x (M. d ) 2 ,p(t,x, z) := p z Jt,x,j). Hence, from ()2.6[) 
V(t, z, y) £ R+* x (R d ) 2 , #(*, z, ») = (*- z, y) = (*, - *„)p(t, 2, y). 

Proposition 31 (Parametrix expansion of the density) With the notations of Sec- 
tion ??, the following representation holds 

oo 

P (t,x,y) = Y,(p® Hir) )( t > x >y)- ( 3 - 4 ) 

r=0 

Proof. Equations (|3.3|) correspond to the forward Kolmogorov equations. Consider now 
the backward equation for p, namely, d s p{s,x,z) = t <l> z p(s,x, z) where t & z stands for 
the adjoint operator of <fr z and the derivatives are taken w.r.t. z. Differentiating under 
the integral we have from (|3.3[) 

rt « r r 

p (s, x, z) p (t — s, z, y) dz 



(P-P)(t,x,y) = J ds^ J 



I ds I [(^$ z p^j (s,x,z)p(t-s,z,y)-p(s,x,z)$ y p(t-s,z,y)]dz=p®H(t,x,y). 
The representation (|3.4[) then follows by simple iteration. □ 

Remark 32 Note that the previous expansion is "formal". The convergence of the 
r.h.s. in (|3.4|) is investigated in the proof of Theorem 3.1 in \KolQ0j and can also be 
derived with the controls of Lemmas \Al\ and \A2\ below. For the sake of completeness, a 
short proof of this convergence is also given in AvvendixWl 
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3.2 Euler scheme 

We consider now, for given JV €E N*, the Euler scheme for equation at the dis- 

cretization times: 

X? = x, X^ +1 = + b(X?Jh + f (Z ti+1 - Z u ) 

recalling h = T/N. 

3.2.1 Proof of Provosition \21\ for the Euler scheme: existence of the density 

For each N £ N*, (X^) i ^^ j^j is a Markov chain. Given the past {X^f = xi, I £ 
[0, i}}, the conditional distribution of the innovations b(X£? )h + f (x^ J ( y Zt i+1 — Zt t ) 
has conditional density ^'(hjO, ■) (with the notation of l|3.1|l . 1)3.21) ). This proves the 
existence of the density for the discretization scheme. 

3.2.2 Parametrix expansion for the Euler scheme 

To give for the Euler scheme an expansion similar to equation ()3.4|) . that will also be the 
starting point for our error expansion, we need to define, for fixed j, k, < j < k < N 
and x, y £ R d additional "frozen" Markov chains (X^)i^y fc j = (Xti,x,y)ley kj- Their 
dynamics is described by 

X# = x, X% +1 = X% + b(y)h + f(y) (Z u+1 -Z u ), ie jj, k - lj. 

Given the past {X^ = x\,l £ [j, i]}, the conditional distribution of the innovations 
b(y)h + f(y)(Zt i+1 —Ztt) nas conditional density p^{h, 0, ■) and, hence, does not depend 
on the past. Note that for the grid points (ti)jg[o n\ the transition densities of the 
solution X s ,x,y of (|3 . 1 f) coincide with the transition densities of the chain X/f XiJ/ for 
N 6 N*, x,y£ R d and s = tj. 

For all < j < k < N, (x,y) € (R d ) 2 , we denote by p N (t k - tj, x,y) and 
P N —~tj, x , y) the transition probability densities between times tj and t k from point 
x to y of the chains X N and X N respectively. In particular, 

p N (tk ~tj,x,y) =p y (t k -tj,x,y) =p(t k -tj,x,y). (3.5) 

Before stating the parametrix expansion ofp N in terms of p N , we need to introduce 
a kernel Hjy that is the "discrete" analogue of H defined in (|2.6|l : 

Hnitk ~ tj,x,y) = |Ljv - L v N ^p N (t k - tj,x,y), (3.6) 

with 

LNVttk ~ tj,x,y) = h l {J P N (h, X, z)(p(t k - t j+ i,z,y)dz - ip(t k - t J + 1 ,x,y)}, 

L v N <p(t k - tj,x,y) = h~~ 1 { j p"(h,x,z)tp(tk - tj +1 ,z,y)dz - ip(t k - t J+ i,x,y)}. 

Note that the previous definitions yield p N (h, x, z) = {^(h, x, z). We also mention that, 
because of the discretisation, there is a slight "shift" in time in the definition of H^. 
Namely we have t k — tj+i instead of the somehow expected t k — tj. 
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Lemma 31 For < j < k < N the following formula holds: 

k—j 

p N (tk -tj,x,y) = ^2(p N ®jv Hfr')(t k - tj,x,y) (3.7) 

r=0 

where in the calculation ofp N ®jy (r-fold convolution) we define 

p (0,x,y) =p (0,x,y) = 5{x-y). 

The proof of this lemma is given in (KMOOj . Lemma 3.6 and does not rely on the 
specific distribution of the innovations. 

( r \ 

Remark 33 With the convention that H N — for r > k — j , equation (|3. T[) also 

writes p N (t^ — tj, x, y) = y^^°_ Q (p N (Sn .ffjy )(tfc — tj, x, y). This expression will often 
be used in the sequel. 



4 Proof of the main results 

In this section, we state in Subsection 14. 1 1 the various points needed to prove Theorem 
1211 The proofs are postponed to Subsection 14.21 As mentioned earlier, the key idea 
consists in comparing the parametrix expansions of the densities p and p N respectively 



given by (|3.4|l and (|3.7[) . In the whole section we suppose that the assumptions of 
TheoremWfXhold. 



4.1 Proof of Theorem I2T1 

For the previously mentioned comparison to be possible we first need to estimate 
a difference between the transition density p(T,x,y) and p d (T,x,y) := ~^^p <S>n 

r>0 

H^ r \T, x,y) which is the analogous of (|3.4|) up to the discrete time convolution (i.e. 
<g> replaced by We refer to l|2.2[). q2.4|). I|2.5|l. Q2.6| l for the definition of operators 

and kernels. 

Lemma 41 (Time discretization) One has: 

(p-p d )(T,x,y) = J2 (7^T)T [P®N (*-**) P d ) ( T ,x,y) 
+h M R Mt i(T,x,y) 

with 

M-l / „ i+1 \ 

J2 (p®n($-$*) P d ) (T,x,y) +\Rm, 1 (T,x, V )\< p a , M {T,y-x). 
1=1 ^ ' 

Then the comparison between p d and p N is controlled with the following 
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Lemma 42 (Comparison of the discrete convolutions) The following expansion 
holds: 



M_1 uk 
(p d -p N )(T, X ,y) = -Y: 



k=l 

M 



(fc + 1)! 



fc+i 



N 



(T,x,y) 



+h lvl R M , 2 (T,x,y) 



where 



R M ,2(T,x,y) = -— I (1-r) 



Ml J 



M 



M+l 



(T,x,y)dr, 



V* e {(*i)ie[i,jv]}, Pr (ti,x,y) = ^2pt<»n hP (U,x,y) ,po =p N 



r=0 



and Vt € [0, 1], pr(t, x, y) = I -p e (rh,x, zffi (t — rh, z,y)dz. 



In particular po(t,x,y) — pp (t, x,y) . Also, 
M-l 



k=l 



fc+1 



P N ) (T,x,y) 



+ \R.M,2(T,x,y)\ < p a ,At(T,y- x). 



Theorem [2T] is then a direct consequence of Lemmas [?T] and [ 



4.2 Proofs of the technical Lemmas 



Proof of Lemma \4l\ We start from the recurrence relation for r £ N* 

p®H [r) -p® N H {r) = \ l {p®H {r - 1) ^j ®H- (p®H( r - 1 A ® N ff] 



+ 



Summing up these terms over r £ N* and using the linearity of d3> and <g>7v we get 
p — p d — p®H — p ®jv + — p d ^j ®n H. An iterative application of this identity 
yields 



(p - /)(T, „) = j[] [p ® - p ® N H] ® N H [r) (T, x, y). 



(4.1) 



By definition, for all fc 6 [1, TV] 



[p® # -p®at H](t k ,x,y) = ^2 I ds f\p(s,x,z)H(t k - s,z,y) 

3=0 ^ ^ 

~P (*j ,x,z)H (t k - tj,z, y)]dz. 



(4.2) 
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A Taylor expansion of the function 0(s,z) := p(s,x, z)H[t k — s,z,y) in the interval 
[tj,s] C [tj,tj + i] gives 



/ 



M-l . ,| . 

e(s,z)-e(tj,z)]dz= V ls jj / a^(r,z)| dz + 



(«-%) 



M r l 



(4.3) 



where Tj (s, 5) = tj +5(s — tj). Note now that —d s p(t — s, x, z) = $p(t — s, x, z), dtp(t — 
s, x, z) — t <l > p{t — s, x, z). Here '<?> = $ z is the adjoint operator of <P where the deriva- 
tives have to be taken w.r.t. z. Hence, $p(t — s,x, z) — t <l>p{t — s,x, z). The same identity 
also holds for p with <P, $ respectively replaced by <P* , We therefore derive 



/ 



d T 0{T,z)\ T=tj dz = d T [p(r,x,z)] \ T=t] H(t k -tj,z,y)dz 

+ I p(-tj,x,z)dr[H(t k -T,z,y)]\ T=tj dz 

= j ' t ®zp{tj,x,z) ($-$*^p(t k -tj,z,y)dz 
- I p{tj,x,z) $*p(t k -tj,z,y)dz 

= I " p{tj,x,z) ($-$*\ p(t k -tj,z,y)dz. 



Iterating the differentiation we get 



d l T e{T,z)\ T=tj dz= j p(tj,x,z) +1 p(t k -t j ,z,y)dz, (4.4) 

where we recall that for two operators A and B we denote by (A — B) k the following 
sum (A - Bf = E) =0 CiA k -i(-By. 

Plugging (|4.3|l and (|4.4|) into (|4.2[) we get 



,/ r_i_i 
[p®H-p® N H](t k ,x,y) = 2^ a + l)! p8lJV (^"^ ) P{ik,x,y) 

1=1 ^ 

+h M R M , 1 (t k ,x,y) (4.5) 

where 

_ I r -i A/ f 1 

R M ,i{t k ,x,y)= {M _ iy Yl I [ h_1 («-*i)J ./ (1 ~ 5)M_1>< 



/ 



J=0 " '3 

3 M 



[p ( r ,x,z)H (t k -T,z,y)] \ T=T .( s S) dsdzd8. (4.6) 
Plugging (|4.5|l and (|4.6|) into (|4.1|l we get 



(p-p d ){T,x,y) = 77-W x £?®JV (*-**) P®n H^ r) {T,x,y) 

1=1 ^ '' r=0 

+h M R MA (T, Xl y) (4.7) 
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with R MA (T,x,y) = ££L (5 M ,1 ® N H^)(T,x,y). 

Now we apply that for a linear operator S and its adjoint S we have p <g) jv <5p = 
t Sp®N p. This gives 



^p®AT p® N H (r) {T,x,y) = 



r=0 



l + l 



P®N^ 
r=0 



P ®AT H 



r) )(T,x,y) 



'+1 



P (T,x,y), 



which plugged into (|4.7[) gives the desired expansion. The stated bound follows by 
application of the estimates given in Lemma [43] below. We only give the proof for the 
first summand, the other terms of the sum over I and the remainder Rm,i(T, x,y) can 
be handled in a similar way. Write 



N-l 

p® N {$-$*fp d {T,x,y) = J2 h 1 P(* 

a n J 



x,z)($-$*y P d {T-t v z,y)dz 



3=0 
Sl+S 2 , 



where in Si (resp. 52) the sum is taken over I\ := {j £ [0, |_ — 2 — Jll ( res P- ^2 — 
{j e [L^J +1,N- 1]}). For Si (resp. 5* 2 ), p d (T - tj,z,y) (resp. p(t } ,x,z)) is non 
singular. From Lemma 1431 below equation (|4.12|) . there exists C := C(T) s.t. for all 
(x,y,z)€(R d f, 

p(s, x, z) < Cj?{s, x, z), p d (s, z, y) < Cp"(s, z, y), Vs e]0, T], 



!(#-*' ) 2 / (T - tj,z,y)\ < Cp y (T — tj,z, y), Vj e h 



p(tj,x,z) < Cp y (t j ,x,z),\/j e h- 



(4.8) 



The semigroup property for p v yields |5i| + IS2I < Cp(T,x,y). One eventually checks 
from Proposition IBll that p(T, x, y) := p~U(T, x, y) < p a m(T, y — x). □ 

Proof of Lemma \42\ Let us denote by ■F^Kz) = / exp(i{z,p))ip(p)dp the Fourier 
transform of a function Introduce now for all u, t, u < t, u,t € {(^')ie[0 N\ }> V £ 

= K l n - L y N )jF(t - u,x, P ) 

p N (h, x, w)p y (t — (u + h),w,p)dw — p y (t — u, x,p). 



Note that in particular according to (|3.6|) . V'(y) — hHpf(t — u, x, y). Taking the charac- 
teristic functions of the densities involved in the above equation, we obtain from l|3.2[l 
and (pTo) that 

F[iP](z) := G,(l) - G,(0) 



with 



G z (t) = exp j(:e,jz} + i(i - u)(B{y),z) + iTh{AB x ' v ,z) 
\{z, s)\ a [(* - u)X f{y) (ds) + ThA\ x > v {ds)} 
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where AB x ' y = B(x) - B(y), A\ x ' y (ds) = \ft x \(ds) - \f^(ds). Note in particular 
that Vr € [0, 1], 



G z {t) = G z (0) x exp rh 



i{AB x ' y ,z) - I \(z,s)\ a A\ x ' v (ds) 
S*-i 



(4.9) 



A Taylor expansion yields T^](z) = h.° { z ] '(0) + ^ /^(l - r) M G { z M+1) (r)dr. 
From 114.911. one derives that for k € N*: 



I G < fc) (0) = ^G,(0) 



:(/Jfl*'V)- / |(^ S }| Q ^A^(d S )l . 



Observe now that G z (0) = J-[6](z), 6(p) := — u,x,p). Using the well-known 
properties of the Fourier transform one gets for all k £ [1, M] 



G ( z k) (0) = h k T 



(*), 



where the operators are applied w.r.t. the x component and the Fourier transform is 
applied w.r.t. the p component of p^(t — u, x,p). Also, in the above writing, we compute 
the Fourier transform for an arbitrary fixed £ £ R d and we then put £ = x. 
Hence, 



m{*) = Efcll i?G{ k H0) + ^ f \l ~ r)^G z M+l \r)dr ■. 

(*) + 



Z^fc=i TT^ 



5 ? - <P y 



TP. 



-T 



Jo a— ) 



A/ 



Af+1 



(*), 



where 



Vr G [0,1], T (p) := / F{Th,x,z)r?{t 



u — rh, z,p)dz. Taking the inverse 



Fourier transform and putting p = y in the above equation, observing that H(t 
u, x, y) — — $ *)p y (t — u, x, y), we obtain 



M-l 



(H N - H)(t-u,x,y) = 2^ (k + i)\ V* ~® J P(* _u ' a: '2/) + 



fc=i 



M\ J 



(1 -r) M (<?>*-<£* ) p T {t-u,x,y)dr. (4.10) 



Af+l 



Recall now that 



oo 

(p rf -p A ')(r, a; ,y) = ^[(?® J v J F/W) - (Wtf^KT,^) 
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where we put yp ® jy H^ N 
identity 



^\ (T,x,y) = for hr > T. Summing over r 6 N in the 



{p® N ff (r) -p®N H^')(T,x,y) = 
((pSjvtf^- 1 ') ®N (H-H N )^ (T,x,y) + 

((p®at # (r_1) ^ r_1) ) ®jv#jv) (T,x,y) 

one gets 

(p d - p N )(T, x, y) = [p d ® N (H - H N ) + [p d - p N ) ® w H N ] (T, x, y). 
By iterative application of the last identity we obtain 

oo 

(p d ~p N )(T,x,y) = Yl {p d ®N (H-H N )j ® N H%\T,x,y). 

r=0 

We get from P~TU)) that for all t £ {(ii)i e [i,j\r]} : 



Af-1 



(p d ® N (H-H N )) (T,x,y) = -J2 ^Y), 



fc=l 



fc+1 



M! 





. /~ ~ \A/+1 


p d ®JV ) Pr 


10 





Eventually. 



ft*" 



fc=i 

AT 



fc+1 



(T,x,y)dr. 



(T,x,y) 



+h M R M ,2(.T,x,y), 



RM,2(T,x,y) = 



(1-r) 



,1/ 



A/+1 



(T,x,y)dT 



Vt € {(*i)i 6 [i,JV]}> pf{t-,x,y) = ^p T ® N (t,x,y) ,po = p N 

r=0 

This proves the expansion part of the Lemma. The bound follows as in the previous 
proof from Lemma [43] □ 

We now state Lemma[43]that allows to control the rests appearing in the expansions 
of Lemmas I41l andl42l Its proof is postponed to appendix [XI 

Lemma 43 Let q > d + 4. For all multi-indices a, b s.t. \a\ + \b\ < q — (d + 4), the 
following inequalities hold: 

I a b d I I a b N I H + I>>l _ 

\DyD x p (t k ,x,y) \ + \DyD x p (t k ,x,y)\ < Ct k a p(t k , x, y), k € [1, nj, (4.11) 
\D%D b x p(t,x,y) \ < Ct-^^pfax^)^ <t<T. 



16 



Also, 3C := C{T) s.t. for all (x,y,z) G (K d ) 3 ,s G]0,T], 

p(s,x,z) < Cp^(s,x,z), p d (s,z,y) < Cp u (s,z,y), 

|(* - $*) V(«, *, »)l < c s -^? y ( s , *, y), 



fe' 

<£-<2> ! 



p(s,x,2) < Cs~ » p"(s,x,z). (4.12) 



5 Extensions and conclusion 

A careful examination of the proofs in the Appendices shows that the absolute conti- 
nuity of A w.r.t. to the Lebesgue measure of S^ -1 can be removed in (A-l) provided 
the function 

£(t,x,y) := J^dJ Rd (J sd l \(P> s )\ ax f(x)( ds ))e*-P (~tj Ifa^r^f^ids) 

x exp(— i{p, x))dp 

has bounded derivatives w.r.t. x up to order q (see Appendix iBl and the statement of 
Theorem 3.1 in [KotOO] ). Also up to a standard perturbative argument, similar controls 
on the density can be obtained when we consider (|1.1[) driven by (Z s + P s ) s >o where 
(Ps)s>o is a compound Poisson process with Levy measure vp(dz) = f(z)dz and 
|/(2)| < 1+ |^jd +)3 ; ft > 0, see Theorem 4.1 in [KolOOj . In that case our main results 
remain valid up to a modification of the remainder. Indeed, it is the smallest exponent 
(or equivalently the largest tail) that leads the asymptotic behavior of p(t, x, y) when 
I a: — y\ is large. Thus p a .Ai(T,y — x) has to be replaced by p m x a l a ,p),M{ r ^i V ~ x ) m 
Theorem l21l Eventually, good controls have been obtained on p for stable-like processes, 
i.e. when the stability index in the generator <Ptp(x) in (|2.2[) can depend on the spatial 
position x, i.e. a turns to a(x) G [a, a] strictly included in (0,2] (see Section 5 in 
KolOO]). Anyhow the processes associated to those generators cannot be approximated 
by a usual Euler scheme and the previous analysis breaks down. The approximation of 
such processes will concern further research. 



A Proof of the controls on the derivatives of the densities (Lemma 

To conclude the proof it remains to prove Lemma l43l The first step is to get bounds on partial 
derivatives of the transition densities p and p. The following estimates generalize the ones 
obtained in [KolOOl . Propositions 2.1-2.3. 

Lemma Al Let q > d + 4. There exists a constant C > 1 such that the following estimates 
hold uniformly for a in any compact subset of the interval (0, 2) and for all < t < T, x,y, z £ 
M d and \a\ < q - (d + 4) 

\D%F(t,x,z)\<J^p*(t,x,z), (A.l) 

\D a z py(t,x,z)\ < — ^ py {t ,x,z). (A.2) 

\z-B(y)t-x\^ 

Remark Al Equation iA.li extends to the stable case what is widely known in the Gaussian 
framework. Namely, each derivation of the density in space remains homogeneous to a stable 
density up to a multiplicative additional singularity of order i" 1 / . 
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Proof. From now on we assume w.l.o.g. that d > 3, the cases d S {1, 2} can be addressed more 
directly. To proceed with the computations, we need to specify a useful change of coordinates. 
Namely, for a given direction f £ R d \{0} introduce for p £ fl d the spherical coordinates 
(p, 1' • ' , <Pd— l)i P = |p| with first coordinate or main axis directed along £, that is 

pi = pcosi?, p2 = P sin i9 cos <p2 j P3 = p sin 1? sin ip2 cos <p3 , . . . 
Pd— 1 = psin i9 sin i/?2... sin 2 cos Vd— 1; 

Pd = psin 1? sin IP2--. sin ipd-2 sm Vd-i i (A. 3) 

i? 6 [0, 7r] , <p; 6 [0, n] , i 6 [2, d — 2], <Pd-i 6 [0, 27r]. Consider then the coordinates (y,T,<j>) 
where r = cos i9 and v = p |C|, with u e R+, r S [—1,1], <f> = (<P2, • • • , <Pd-i) S [0, ir] d ~ 3 X 
[0, 27r]. In the following we write p = piy, r, <j>) for the previous r.h.s. in I IA.31 written in these 
new coordinates that is 



PI 


= ICI 


VT, 


p 2 = icrvi-r 2 ) 1 / 2 


cos 1^2 , 


P3 


= ICI 


_1 v(l 


— t 2 ) 1//2 sini^2 cos<p3, . 




Pd-1 


= ICI 




— r 2 ) 1//2 sin^2--- sinipd- 


-2 cos <p_-_ 


Prf 


= ICI 




— r 2 ) 1//2 siny>2... sinipd- 


_ 2 sinipd-i 



(A.4) 

and p(r, <j>) = p(|C|, r, <f>). 

Without loss of generality we suppose B(y) = 0. The first step consists in differentiating 
w.r.t z the inverse Fourier transform for p y (t,x,z) 

p"(t,x,z) = f exp(-t f \{p,s)\ a X f{y) (ds)\cxp(~i{p,z-x))dp. 

(A.5) 

For z = x, 112.31 1 and standard computations directly give estimate iA.lll . Thus, in the following 
we also assume z/i and use the previous spherical coordinates (v, t, <j>) derived from J A .3I > 
setting £ = z — x as the main axis. We obtain: 

D?pU(t,x,z) = r-r-7 / dvv^+^x 

Z V ' {2TT) d \ Z - X\\ a \+ d Jo 

f At I d(f>9(v,r,\a\)exp\-t V f \(p,s)\ a \ f(y) (ds)\x 

J-l /[0,7r]<J-3x[0,27r] I \* ~ X\ is- 1 - 1 J 

\a\-a 1 +d-3 

T ai (l-T 2 ) 2 h(<p,a) , (A.6) 
where p = p/ \p\, a = (ai , ad) S N d and 

${v,T,\a\) = (-l)l a l/ 2 cos( V r)I| a | even + (-l)Cl a l+ 1 )/2 S m{vr)I M odd , 

h((p, a) = {(cos <P2)° 2 (sin ip2 cos ip3 ) a3 X ... X (sin <p_... sinipd-2 cos tpd-l) a<i ~ 1 
X (sinip2---sinipd-2sin¥'d-i) ati } X V(<f>), 
V{4>) = (sin(p 2 ) d ~ 3 (sinip3) d ~ 4 X ... X (sin ipd-3) 2 sin Vd-2- 

We consider, first the case \z — x\ ft 1 / 01 < C, for a sufficiently small positive constant C. In 
this case we expand the trigonometric function ^(v, r, \ a\) in JA.6I in power series and change 

the variable of integration to w in each term. This gives for all tgN, 



\i\ + d 1 dm 4-1 < i jjV m \ t 1 / c 

tS— L m=0 \ Zm + Ha\ odd)- \ £ 

1 (-_-n(M+I|a| odd)/2 

jM I r>. . — > _____ 
(2n) d 



. C H = , n _ w , (A.7) 
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where Vm e [1, kj, 

e™ = / dw dr I d<f>cxp \ -w a / \(p, s)\ a \ f(v] (ds) } 

JO J -I y[0,7r]<i-3x[0,27r] I J S^ 1 J 

w \ a \ +2m +d-l M even x r «l+2m+I| | odd ( X _ T 2 ) ~l + d ~ 3 a ), 



k&jl ^l^-^iyC^+llal odd 

(A.8) 



" (2(fc + l)+I, , odd )I I tV- 



To simplify the notations we omit the dependence of the coefficients of our expansions on the 
direction £ = z — x. From (A-l), (A-2) and 1 12,31 1 one then derives the following bound: 

\a\ + 2m + d + I| a | odd 



j?( m+ ai + 1+I|a| ° dd , |a| - ai+d - r ). (A.9) 



Here denotes the area of the unit sphere S d ~ 2 and B is the /3-function. Note that the 

modulus of each term in the expansion iA.7\ serves as an estimate of the remainder in a finite 
Taylor expansion. From IA.7\ we have 

w.-.*) = (4 a| (^?) lModd + M") ■ (A.10) 

Recall that we are considering the case < C. By Proposition 3.1 (i) from IKolOOl for 

some C depending on C, C~ x t- d l a < p"(t,x,g) < Ct- d / a . Hence, equations ifATOl . fOt . 
fOl yield 

\D a z P>(t, x,z)\< -^-rpUft, x, z) < CC , . p»(t, x, z). (A. 11) 

~ t ^ |z-z|H 

To estimate D°p«(t,x,z) for |z - x| /i 1 /" > we proceed as in Proposition 2.3 of 

KolOO]. This gives the following representation D^p v (t, x, z) = [D^p v (t, x, z)] 1 + [DJp ! '(t, x, z)] 2 
with 

i poo rl 

[Dlp»{t,x, Z )\i = _ jf < ^pW+«'- i y_ i dT*(p|*-«| ) T ) |a|) x 

/j(r) / #cxp{-tp a 9A (r,</>,j/)U(0,a), j = 1,2, (A.12) 

g x Ar,4>,y) := / |(p(r, 0), s>| Q A /(y) (ds), 
using the notations introduced after I IA.31 . Here 

|q|- ai +d-3 |a|- ai +d-3 

/l(r) = r a ± (1-t 2 ) ^ XW, / 2 (r) =r a i (1-r 2 ) = (1 - X (r)) 

where x( r ) is a C 00 even truncation function M — >■ [0, 1] that equals 1 for \t\ < 1— 2e, and for 
\t\ > 1 — e for some £ £ (0, j). Because of the symmetry in r, it is easy to see that the integral 
in is non-zero only if ai and \a\ are both even or odd. Expanding the exponential at 

order 2 in HA. 121 1 and making the change of variables p\z — x\ = v we get 

[DSjFfcx,*)]! = C| °' +d £ A 6 - 1 (r^T' ( A - 13) 
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where C\ a \ is defined in (I A. 7l l and for m £ [0, 1] 



and 



<fim (t) 



(-l) m / F;^{v)v^+ ma+d - 1 dv, 



a even 



Rc- 



Jm] 



r 



exp(— ivr)ip m (r)dr 



Mr) [ gFAT.hvMhaW. 

V[0,vr] d -3x[0,27r] 1 



./o Jo 



F 2s(. v ) = evenRe-I| a | odd Im] 



r 



exp(— ivr)tp2 s{t)<1t 



<P2,s(t) = /i(t) 



/ 



! x[0,2tt] 



9\ f (?) 0: V) ex P 



-5t 



\z - x\ 



9\f(.T,(f>,y)\ 



xh(<f>, a)d<j>. 



To extend the integration to R in the definition of (-Fm' ( u ))me[0,11 1 ^"J^)' we s i m Ply use 
that the functions (Vm) m g[o,i]i V2,<5 have compact support in r. However, to check that 
the coefficients (&m)n»e [0,2] are well defined, we have to equilibrate at infinity the term in 
^ v \a\+ma+d— 1 ) m g|Q j2 ) . This can be done computing iterated integration by parts in r in 

the definition of (F„ ' (v )) m g[o,i] i F^l^)' Namely, if ip m {r),m = 0,1, and ip2,s( T ) are C 9 
functions of r with compact support and q > |a| +i + d > |a| +2« + d, performing <jr integrations 

by parts w.r.t. r one derives that the coefficients (&m') m g[o,2| are weu defined. Let us now 
check that assumption (A-l) implies that ipm(r), m = 0, 1, and ip2 $(t) are C q functions of r 
with compact support. Indeed, for the unit vectors p(r + At, <j>) and p(r, <f>), from elementary 
algebra there exists an orthogonal matrix A := A(At) s.t. p(r + At, <j>) = j4p(t, </>). Hence, if 
^f(x){ds) = s)ds where has the previous smoothness one can show 



lim 

Z}t->0 



lim 

/4t->0 



SAy (l" + ^T, <£, - 9Ay (T, 0, z) 

At 

f sd - 1 {\(p(T,<f>),A*s)r - \(p(r,4>),s)\ a }X f(x) (ds) 



At 



[ |(p(r,0), S >r lim 



[e(i,As)-e(x,s)] 



(/,S 



where /3(t, </>, s) is C°° function in t uniformly bounded in (r, </>, s) in our region. The process 
can then be iterated other q — 1 times. 

Thus all coefficients (&m') m g[o,2] are weu defined. 

Next, analogously to Proposition 2.3 in IKolOOl (where the case |o| = was considered) 
and with the same rotations of the integration contours for a £ (0, 1], a S (1,2), we obtain 
for all k 6 N 



[D a z P(t,x,z)] 2 



C 



\z - x\ 



\a\+d 



k 

V — e H 



m = 



\*-x\ a 



n 2,k+l 



(A. 14) 



2[I|a| cvonRe-I|a| oddMl/ dr l dcf>h((j>, a) (-g\ . (t, 4>)Y 

Jl-2e J[0,7r] d -3x[0,27r] 



X exp(- 
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and l-R^I+ll — (fc+j 1 )! ( | x — z | " ) ' Note that the coefficients cj^' are also well defined 
because r does not approach zero (recall that 1 — x( T ) 7^ |r| > 1 — 2e). Precisely 
|cm'| < 2A d _ 2 C!p(l - 2£)- am + d +l a lr (am + d+ |g|). 

Now the sum of expansions J A .131 ) and <A.14I gives the expansion for D^p y (t, x, z). Note 

that by construction, the first coefficient + c|, a ' does not depend on the spectral measure 
and it vanishes when the spectral measure is uniform (that is C\ = Ci = 1 in 112, 311 ). 
This can be shown by means of representations involving Bessel and Whittaker functions and 
the same rotations of the integration contours as in Proposition 2.2 of IKolOOI . see Appendix 
[C]for details. Thus, for all k S N*, we get a representation 

^^-j7^AtM aj (^r +<i }' (a - i5) 

where di^j = &m + 4^' with b m = for m > 3 and l-R^J < + ' Now ' from 

Proposition 3.1 (ii) in IKolOOI > 0. Equation flA.151 yields 



|z-x| d V|2-^r / " 2\z-x\ d+a 

for sufficiently small C. Hence, we have 

\DlpK{t,x,e)\ < CC]a \ , < 2_^p»(t )X)Z ) 

|z-x|l a l |z-z| d + Q | 2 -x|l a t 

< fSrp"(*.*.*). (A.i6) 

recalling that \ z _ x \ < C f° r the l as t inequality. W.l.o.g. we can assume G < 1. It remains to 

consider the case |x - ^l/i 1 / eJC.C^ 1 ^ /(C). It follows from fOt that |jg - z| d p»(t, x, *) 
and \z — x\d+\ a \D®'p y (t,x,z) are continuous functions of |x — z\/t 1 / OL . Since the stable density is 
also strictly positive, we deduce that there exists C s.t. on 1(C), \D'^p y (t, x,z)\ < ^\ a \+d ^ 

1 — C || a | P^fti x i ■ 2 ) ^ ^~M/t» P^ft' z ) which concludes the proof. □ 

Lemma A2 Let <jr > d + 4. There exists a constant C > 1 s.t. the following estimates hold 
uniformly for a in any compact subset of the interval (0, 2) and for all < t < T, x, y, v £ R d 
and \a\ +\b\ < q- (d + 4): 

\DlDlH(t,x,y)\ < -^pft.x.y) (l + """^ = *' ) , (A.17) 
| x, a: + < Cp(t, x, x + v) fl+ mm ^' l v l j , (A. 18) 



D;r>*p(t,*,»)| < , „,_ A " .ii.um K'.'.y)' (A.19) 



|2/- B(yjT- x\M+\ b \ 
Proof. Inequalities J A. 1 71 ) and flA.181 follow from the representation 



H{t,x,y) = (B{x)-B{y),V x p{t,x,y)) + -±- ( \p\ a ( \(p,s)\ 

{\ f(y) {ds) - \ f(x) (ds)) cxp \-t\p\ a / \(p,s)\ a X f{y) (ds) 
I Js d ~ 1 
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exp{-i(p, y - B{y)t - x)} dp (A.20) 

analogously to the proof of Proposition 2.3 in IKolOOl , see also Appendix FBI where J A . 171 ) is 
proved for \a\ = \b\ = 0. Inequality (I A . 18ft contains in (3.23') p. 748 of that reference. Inequality 
< l A .191) can be derived following the proof of Lemma IA1I □ 

The proof of Lemma 1431 can then be achieved from Lemmas IA1I and IA2I adapting the 
arguments in Appendix[B]concerning the control in terms of the frozen density for the "formal" 
series appearing in 13.41 . See also the proof of Theorem 2.3 in IKM02I or Theorem 3.1 in [KolOO . 



B Control of the parametrix series of the density 

For the sake of completeness we provide in this section a complete proof of the control for the 
r.h.s of 113.411 under our standing Assumptions (A-l)-(A-3). 

We first sum up in Proposition IBll the various estimates needed to control the convergence 
of 1 13.41 1 following the proof of Theorem 3.1 in IKolOOl . namely Proposition 3.1 and its corollary, 
Lemma 3.1 and Propositions 3.2-3.3 of that reference. These estimates can also be directly 
derived from the computations of Appendix ICl 

Proposition Bl For all K sufficiently large, there exists C > s.t. the following estimates 
hold uniformly for a in any compact subset of (0, 2), for all x,y, z £ Mr and for all t S (0, T], 

c -i t -d/ a < x> z ) < ct~ d/a , \x-z\< Kt 1 / 01 , 



\x-z\ d + a \x-z\ d + a ' 
p z (t,x,z) < Cp"(t,x,z). 

Also, there exists C > s.t. V(t,x,y) £ [0, T] X '~® d ^ 2 



dz min(l, \z\)p"(t, 0, z) < Cf , uj : = min(l, 1/a). 
For all s e (0, t) 

j dzp z (t — s, x, z) min(l, \y — z\)s~ 1 p u (s, z, y) 
< C^- 1 min(l, \y - x\) + s" -1 )^^, x, y), 
fdzmin(l,\z-x\ W (t-s,x,z) 
X min(l, \y - z\)s~ l p"{s, z, y) < OV" V(*> x, y), 
J dzp z (t — s, x, z)p y (s, z, y) 
< Cp»(t,x,y). 

Introduce now for a given bounded measure r) on S d_1 the function 

<*»(*.«.*«,)> ■■= ^L d ^f sd j^ 

xcxp (~t\p\ a J^ ^fas^Xf^ids^j exp(-i(p,*». 

With this notation and HA. 201 ) we get 

H(t, x, y) = (B(x) - B(y), V x p(t, x, y)) + (<Px f{y) ~ VA /(<r) )(*, y-x- B(y)t, X f(y) ). 

Under our standing assumptions, the mean value theorem yields IvXfjy-j ~ fkf;^ )(*> V ~ x ~ 
B(y)t, Xf( v ))\ < min(l, \y — x\)ip Vxy (t,y — x — B(y)t, A^( y )), where rj x , y is a bounded measure. 
Now Proposition 2.5 in IKolOOl states that for a bounded measure n, 

tp v (t,z,\ f(y) ) < Ct-V(*,0,2). 
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From Lemmas IA1I lA"2l and the above controls one deduces \H(t, x,y)\ < Cp 11 (t,x,y)(l + 
f. -1 min(l, \x — y\j) := Cv(t, x, y) (which actually gives HA. 1711 for a = b = 0). 

Introduce now /3 o ip{t, s, x, y) = / /3(t — s,x, z)ip{s, z, y)dz, i.e. o is the spatial part of the 

convolution operator ®, and set v(t,x,y) := tv(t, x,y). From Propositon IB1I one derives 

po v(t, s, x, y) < C(v(t, x, y) + p(t, x, y)) 
v o v{t, s, x, y) < C(v(t, x, y) + + (t — s)^)p(t, x, y)). 

Recalling \H(t,x,y)\ < Cv(t, x,y), integrating the above inequalities one gets: 

|p® H\(t, x,y) < C{v(t,x,y) + t^p(t,x,y)), \p®H® H\(t,x,y) < C 2 ^ (p + v)(t, x,y). 

An induction yields, for all k S N*: 

\p^H^ k 1\(t,x,y) < — (p + v)(t,x,y), 
{k\) 2 

\p^H( 2k+1 1\(t,x,y) < — -(tv(t,x,y) + t"p(t,x,y)I a>1 ), 

fc!(fc + I Qe(0 ,i])! 

and the the required control, i.e. p(t, x, y) < Cp(t, x, y). The controls on the derivatives can be 
proved in a similar way, up to suitable rearrangements of the variable of integration, see p. 747 
and 748 in IKolOOI . Also the whole proof can be carried out for p d ,p N . □ 

Remark Bl To conclude, note that by arguments similar to those used to prove Proposition 
\B1\ one gets 

\H® H(t,x,y)\ < Ct u '- 1 p(t,x,y), 
which turns to be a useful estimate to derive 14.81 following the above proof. 



C Additional computations concerning the derivatives of the density 

In this section we prove that b^' + c{, a = 0, justifying that the first index in 1A.15I is o 

Odd dimensions d 

From the definitions in the proof of Lemma IA1I . it is enough to show 



I|o| ovcnR e - I|a| odd lm | 



■J^ + ^dv 



+2 



If f cxp (-mit) fi{r)di 
Uo Uv. 

f 1 (ir)-( d +l Q l)r(d+|a|)/ 2 (r)dr) = 0. (C.l) 



Note that since d is odd, if \a\ is odd j-( d +H) = (-1) 2 and if |a| is even i ( d +H) 
writes 



d-l+|a| 

i~ 1 (— 1) 2 . Hence, the contribution of the second term in lC.lt vanishes and the condition 



[l|o| eve n Re-I| a | oddH [J ^exp(-ror)/i(r)dr v^ +d 1 dv\ 



(C.2) 



Denote Gi(v) = J m exp (— ivr) fi(r)dT. Remind that \a\ and a\ have the same parity, see proof 
of Lemma IA1I 
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a) For even |a|,ai, Gi is even and belongs to a Schwartz space of functions. Since d is odd, by 
inverse Fourier transform, Equation JC 2t reduces to 



dv 



rOG i p 

/ Gi(y)vW +i - 1 dv = - / G 1 {v)v^ +d - 1 
Jo 2 y H 



-/{'"'^''(O) =0. 

The equality y(' a ' +d 1 '(0) = follows from the Leibniz differentiation rule for the product 

H-oj+d-3 

t" 1 X (1 - r 2 ) 2 and the identity |a| + d— 1 = ai + (|o| - ai + d - 3) + 2]. Thus 

IC.2I holds in this case. 

b) Analogously, for odd |a|, a\ , -Im(Gi (i>)) is odd and belongs to a Schwartz space of functions. 
The function (-ImGi(i>)i>l a l+ d ~ 1 ) is even. Thus 



f 
Jo 



(-lmG 1 (v))v^+ d - 1 dv= ( i) ^' +d / 1 (|a|+d ~ 1) (0) = 



10 2 
for the same previous reasons and equation llC.21 holds in this case as well 



Even dimensions d 



We assume in this section that the spectral measure is uniform with Ci = G2 = 1 in (A-l). 
For \a\ and ai = 2m even, equation JA.6I 1 can be rewritten as 



(_l)kl/2^a m 



(2ir) d \z - x|M 



(C.3) 



3=0 



X J™ dwl a l+ d - 1 exp (-t \ J 1 (1 - t 2 ) n i - 1 ' 2 cos(vr)dr 



where A°; 



~ J[0,7T]<*-3x[0,2,r] h ^ a W and N 3 



| a |_ ai + d _2+2j [ j g |0, m]. Now recalling 



the definition of the Bessel function J n {z) := 1 _ Ji(l — t 2 ) n 1//2 cos(zt)dt which 

well defined for n > 1/2 on C\(— 00, 0), we get 

(_l)|a|/2^4a m , 



(2-K) d \z - x\\ a \+ a 



3=0 



./() 



J | + q 1 +d-2j / „" \ 

x/ v 2 exp ^-t . j (u)<fo 

f 1 \|a|/2 Aa m , 

I l ) A d~2 V^< 1. 

3=0 



d-2 y-(„ 1) 2m-3 C ^ r(7V + 1)2^+1/2 (C.4) 

(27r) d |^-x|i [l i+ d +-r 3 r 



xRe 



r exp (- 



\z-x\ a 



exp 



2 J 4 



where W ,„(z) = c * p . ( ~ z /f / °°[t(l + t/.z)] n ~ 1 '' 2 e- t dt, n > 1/2, 2 6 C\(-oo,0), is the Whit- 
taker function and for z > 0, 



J n {z) = 2Re 



7 =exp(i(n + -)7rt)Wo, n (2tz) 
/27TZ 2 2 
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(relation (2.10) from IKolOOI I and 7VJ = |a|+ " 1+ 2 ' i ~ 2j ~ 1 , j 6 [0,mj. For a G (0,1], from 



Cauchy's theorem we can change the integration path in IIC.4I to the negative imaginary half 
line. Setting then v = —iS; we obtain 

(-l)\a\/2y^a m 

D a 2 p»(t,x,z) = ——f f=g-i y)(-i) am - j c&r(jv j - + i)2 Ar J + 1 /2 

(2n) d \z - x\^ +d f^ 2 



Jo 



-t I exp ™ ) ) W 0tNj (2S)t N 'i d£ 



Recalling the definition of Wn.iVj > we conclude expanding the exponential in power series that 
the first term is 0. 

For a £ (1,2), using the same arguments we can rotate the initial contour through the 
angle — 7t/(2q). Setting then rj = e l 2&v we get 

f_l~)|f|/2^a m 

D a z p«(t,x,z) = ——f f=g-i yV-i) am - J c£ l r(iv j + i)2^ +1 /2 

(2n) d \z - x\^ +d f^ 2 



xRe / exp I it- — - — ■ h ( -Nj + - ) m - — (N 1 , + 1) 



Jo 



\z-x\ a \2 



2a 



XW , N] (2„ exp { J. )r? ^ drj . 



2a 



Taylor's formula for cxp(tf \ z ^ x \ a ) yields for the first term, Vj G [0, m], 

'i7r(a- 



I£ := Re < exp 



Ijvi + ijwi-— (JVi + i) 

2 J 4/ 2a J ' 



/■ OO 

J W , Nj (2r 1 ( 



cxp^^_^)^ dr?} 



7T (a: — 1) 



2 a 



At last, we rotate the contour through the angle 

obtain II =Re{-i(-iy- m / oo W 0l JV i (20e JV ^} = 0. 
For \a\ and ai = 2m + 1 odd, 



Setting £ = f?exp 



in (ol — 1) 



(_l)(|a| + l)/2 J 4a m 

DW(t,x,») = 7T-rnrY(~ 1 f m ~ j cL 

z V ' (2vr) d 2 - z M+ d f-! v ' m 
X (W^" 1 exp ^ (1 - r 2 )^- 1 /^ sm(vr)dr 



(_1)(I»I+1)/2A^_. 
(27r) d |z - ie|M + 



_^2 ^ ( _ 1)2m+1 _, c4 / drala | +d -2 cxp ( t 



|z - x\ c 



x y 1 (1 - r 2 ) JV 3- 1 / 2 Td(cos(r;T)) 

(•_l~)(M+l)/2y|a "i r oo / ,,a \ 



(27r) d |2 - z| 

X ^ cos(„t)x [(1- 2iVj)r 2 (l -r 2 ) JV J " 3 / 2 + (l-r 2 ) JV ^ 1 '' 2 ] dr. 



The above integrals have the same form as in JC 3t and can be estimated similarly. 
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